Tunable Kondo physics in a carbon nanotube double quantum dot 
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We investigate a tunable two-impurity Kondo system in a strongly correlated carbon nanotube 
double quantum dot, accessing the full range of charge regimes. In the regime where both dots con- 
tain an unpaired electron, the system approaches the two-impurity Kondo model. At zero magnetic 
field the interdot coupling disrupts the Kondo physics and a local singlet state arises, but we are 
able to tune the crossover to a Kondo screened phase by application of a magnetic field. All results 
show good agreement with a numerical renormalization group study of the device. 
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A quantum dot coupled to two leads can be consid- 
ered as an experimental realization of the single-impurity 
Anderson model [l|. When the dot contains a single un- 
paired electron, the Anderson model accurately describes 
how, below a characteristic temperature Tk, correlated 
electron tunnelling between the quantum dot and the 
leads results in coherent screening of the electron spin 
• The combined system of electrons on the quantum 
dot and leads forms a spin singlet, a phenomenon known 
as the Kondo effect |5|. Likewise, two tunnel-coupled 
quantum dots should amount to an experimental realiza- 
tion of the two- impurity Anderson model Q. Here the 
physics is much richer, particularly in the regime where 
each dot contains an unpaired electron. In this case, 
a competition now arises between the tendency of the 
conduction electrons on the leads to screen the spins on 
the quantum dots, and the antiferromagnetic exchange 
coupling J between the two localized spins. The former 
favors formation of a Kondo singlet between each lead 
and the dot to which it is coupled, while the latter fa- 
vors a local singlet state. The resulting groundstate of 
the system depends sensitively on the relative strength 
of the interactions, an understanding of which is impor- 
tant and believed to underlie the electronic properties of 
a wide range of strongly correlated materials, including 
spin glasses and heavy fermion compounds Q . 

The essence of this competition is captured by the two- 
impurity Kondo model |7H10|, which describes the low- 
energy physics of the two-impurity Anderson model in 
the absence of charge transfer between the leads, and 
famously contains a quantum phase transition at the 
boundary of the local and Kondo singlet phases where 
J ~ Tk- While this has attracted considerable exper- 
imental attention 11 -13|, observation of the transition 
has remained elusive. This is perhaps unsurprising as 
charge transfer between the leads, absent in the two- 
impurity Kondo model, is necessarily present in experi- 
ment if a conductance is measured. As is well known the- 
oretically, this transforms the quantum phase transition 
into a crossover, such that the ground state of the system 



is always a Fermi liquid [ToL 14 1; although remnants of 
the transition are evident in a strong enhancement of the 
zero-bias conductance in the vicinity of J ~ Tk flET- Y!\. 

An understanding of the transport properties of a re- 
alistic two-impurity system, such as a double quantum 
dot (DQD), thus requires that charge transfer between 
the leads is taken into account. This is achieved in the 
present work, where we present a study of a tunable car- 
bon nanotube DQD in the strongly correlated regime, 
and use a numerical renormalization group (NRG) study 
of the two-impurity Anderson model to describe the de- 
vice. We show that in the charge regime where both dots 
have an unpaired electron, the ground state of the device 
is a local singlet phase with suppressed Kondo correla- 
tions. The ability to tune the exchange coupling and 
the Kondo scales allows one, in principle, to crossover 
from the local singlet to the Kondo screened phase. In 
our device the onset of charge fluctuations prevents the 
crossover being seen cleanly at zero magnetic field. We 
have, however, been able to observe it at finite magnetic 
field, consistent with recent theoretical predictions [3 ■ 

The device we consider is a single- walled carbon nano- 
tube on a degenerately doped Si/SiO? substrate con- 
tacted by Au contacts, see Fig. 1(a) [18J. A central gate 
is used to introduce a tunable tunnel barrier, separating 
the nanotube into two quantum dots, which can be in- 
dividually addressed by two additional side gates. The 
stability diagram of the device is shown in Fig. 1(b). The 
effective electron number of each charge regime is indi- 
cated by the ordered pairs (n, m). The large-small-large 
alternation in the stability diagram reflects the two-fold 
spin degeneracy of the device, and allows us to establish 
unambiguously the parity of the electron number (even or 
odd). The orbital degeneracy of the nanotubes is broken, 
most likely as a result of K — K' mixing 0; 0] • 

From an analysis of the stability diagram we are able 
to extract the on-site charging energies U ^2.5 meV 
for both dots, the electrostatic coupling energy U' ~ 0.6 
meV and the interdot tunnel coupling t ~ 0.4 meV. The 
widths of the Coulomb blockade peaks in the stability 
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diagram allow an estimate of the coupling T between the 
dots and the leads. In the non-interacting limit the full- 
width at half maximum is equal to 2r, but the peaks 
are broadened by a further factor of ~ 2 by many-body 
scattering processes 2l|. Bearing this in mind, we ob- 
tain T R ~ 0.23 meV and T L ~ 0.12 meV for the right 
and left dot respectively. For these coupling strengths 
U /Y v ~ 10, so we might expect co-tunneling processes, 
and thus the Kondo effect, to be experimentally observ- 
able. To investigate the presence of Kondo correlations, 
we measured the differential conductance in all charge 
regimes, see Fig. 1(c) and 1(d). When both dots con- 
tain an even number of electrons, the Kondo effect is in- 
operative and conductance is suppressed and featureless 
around source-drain bias V s d = 0. On the other hand, 
when one of the quantum dots contains an odd number 
of electrons, we observe a pronounced zero-bias conduc- 
tance peak. The evolution of the differential conductance 
as electrons are added to the right dot, keeping the ef- 
fective electron number of the left dot fixed at N = 2, is 
shown in the top-most panel of Fig. 1(d). The charac- 
teristic appearance of a zero-bias peak when the electron 
number is odd is a clear indication of Kondo physics in 
the (0,1), (1,0), (1,2) and (2,1) charge regimes 0l3|, the 
spin of the singly-occupied dot being effectively Kondo 
screened by the lead to which it is coupled. 

The behavior, however, is markedly different in the 
center of the (1,1) charge regime where the electron num- 
ber is odd for both dots. As shown in Fig. 1(c), the 
zero-bias conductance is suppressed and a double peak 

mm. 



The behavior we 



structure arises at finite bias 
observe is characteristic of a strongly correlated DQD 
when J > T£,T£ where T£,T$ denote the Kondo 
scales for the left and right dots respectively. These 
energy scales are readily estimated using the parame- 
ters obtained above. In the center of the (1,1) charge 
regime the exchange coupling J ~ 4i 2 /{U — U'), which 
yields J ~ 0.34 meV. The Kondo scales can be estimated 
roughly using the Haldane expression [22| 



T V K ~ Vr7f7exp[^(e, + U)/2Y V U] 



(1) 



where e v (v = L, R) is the level energy of dot v relative to 
the zero-bias Fermi level of the leads. This yields T§ ~ 
10" 2 meV and T£ - 10~ 4 meV in the middle of the (1,1) 
charge regime, such that J 3> T^,T^?. 

The above interpretation of the measurements is 
strengthened by an NRG study of our double quantum 
dot. The calculated stability diagram using the exper- 
imental parameters, see Fig. 2(a), reproduces all key 
features of the experiment and allows for a comparison 
between experiment and theory (for further details see 
supplementary material [18j). Finite-bias conductances, 
NRG calculation of which is inevitably approximate (see 
[IH), are also in good agreement. In particular the dou- 
ble peak structure in the (1,1) charge regime is repro- 
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Figure 1: (color online) (a) Schematic of the carbon nanotube 
DQD device, (b) Stability diagram of the DQD measured at 
T ~ 60 mK. The Si backgate voltage V bg = -1.7 V. The 
barrier gate voltage Vbar = mV. (c) Differential conductance 
in the centers of the various charge regimes of panel (b) as 
indicated by (n,m). The curves in the rightmost panel are 
offset in steps of 0.03 e 2 jh. (d) Differential conductance along 
the lines indicated in panel (a). 
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Figure 2: (color online) (a) NRG calculation of the stability 
diagram, (b) Differential conductance in the centers of the 
various charge regimes of panel (a) as indicated by (n, m). (c) 
Differential conductance along the lines indicated in panel (a). 
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duced in the calculations, see Fig. 2(b). It can be under- 
stood physically as a non-equilibrium Kondo effect (23|. 
While the formation of a local singlet state suppresses 
Kondo correlations at zero-bias, interlead spin-flip tun- 
neling becomes possible and Kondo correlations are par- 
tially restored when V s d is comparable to the exchange 
energy J separating the atomic-limit singlet ground state 
from its triplet excited states. As demonstrated below, 
a unique feature of double quantum dots is that Tj^,T^ 
and J are all tunable by varying the dot level energies. 

We discuss first the behavior along diagonal 2 in 
Fig. 1(b), for which the differential conductance is shown 
in Fig. 1(d), middle panel. Along this diagonal, the dot 
level energies cl,r are decreased while their difference 
(or detuning) ex — en remains zero. This allows us to 
tune the Kondo scales of the two quantum dots, as these 
strongly depend on the dot level energies, see Eq.([T]). In 
the center of the (1,1) charge regime, the Kondo scales for 
both dots are at a minimum. Moving away from the cen- 
ter, the Kondo scales thus increase. However, before we 
can access the Kondo screened phase we reach the edge 
of the (1,1) charge regime, moving into a mixed- valence 
regime. Here charge fluctuations are significant, and a 
two-impurity Kondo model description is no longer valid. 
The inability to access the Kondo screened phase this way 
can be understood given that for our device J ~ ^l,r, 
so to observe the crossover requires T^' ~ Fl,r, which 
occurs only when ex.x?| < T^^, i.e. in the mixed- valence 
regime. A lower J is therefore required in order to ob- 
serve the crossover cleanly in this way. While we were 
able to control the tunnel coupling t experimentally (see 
[IH), sufficiently small values of it were not reached here. 

We now focus on the behavior along diagonal 3 in 
Fig. 1 (b) for which the differential conductance is shown 
in Fig. 1(d), bottom panel. Along this diagonal the de- 
tuning (ex — en) is increased while ex + ex? is constant. 
As indicated by the dashed white line in Fig. 1(d) the 
peaks observed in the differential conductance move fur- 
ther apart as the magnitude of the detuning is increased 
and the conductance becomes highly asymmetric in bias. 
This increase in peak splitting is due to J becoming larger 
with positive detuning as the (1,1) sin glet state becomes 
closer in energy to the (0,2) singlet [24j. The ability 
to tune the exchange energy by varying the detuning 
is essential in spin-based quantum information process- 
ing schemes using quantum dots (25| . Importantly, the 
present data shows that this tunability can also be used 
as a probe of Kondo physics. 

The consequences of varying the exchange J are in- 
vestigated further in Fig. 3 by application of a magnetic 
field (B) perpendicular to the nanotube axis, for vari- 
ous values of the detuning along diagonal 3 in the (1,1) 
regime. At finite B the observed zero-field peaks split 
into three components, Fig. 3(a) (of which the innermost 
peaks are most easily resolved, see Figs. 3(a,b)), consis- 
tent with a local singlet ground state and triplet excited 
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Figure 3: (color online) (a) Magnetic field dependence of the 
differential conductance in the (1,1) charge regime for various 
values of the detuning as indicated in Fig. 1(b) (black dots 
along line 3), with the left plot for the center of the (1,1) 
regime, (b) NRG results corresponding to the left and right 
panels in (a). Dashed lines show the threefold splitting of one 
zero-field peak, (c) Schematic of crossover between LS and 
KS phases, see text, either by varying B (solid arrowed line) 
or J (dashed arrowed line), (d) Zero-bias conductance vs B 
for different detuning, along dashed white lines in panel (a); 
both measured (left) and calculated (right). Individual curves 
from left to right correspond to (a) from left to right. 



state, separated by ~ J at B = 0. As B is increased, 
the energy difference (~ J — gunB) between the singlet 
(S) and the lowest-energy triplet state, TL, decreases. 
At a field B c ~ J/gfin these states are near-degenerate. 
A strong zero-bias conductance peak is then observed, 
due to Kondo screening of the S-T_ pseudospin system 
(26M28| . here by two channels. [In our tunnel-coupled 
double QD both the S and T_ states arise from the (1,1) 
charge configuration, in contrast to S-T- crossings ob- 
served 0, 0, [29|, [30I in single QDs with two 'active' 
levels, where the relevant singlet is a configuration with 
the lower dot level doubly occupied [3l|, and conse quen t 
Kondo screening arises by a single effective channel [26|.] 

From the perspective of the two-impurity Anderson 
model, and its experimental realization in our device, 
this conductance peak is the signature of the finite-field 
crossover from a local singlet (LS) phase for B < B c , to 
a polarized Kondo screened (KS) phase for B > B c : and 
amounts to a finite- B continuation of the zero-field con- 
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ductance peak at J ~ Tk (l4 - 17|. For a L-R symmetric 
two-impurity Kondo model, it was recently shown (l4| 
that the well known zero-field quantum phase transition 
between LS and KS phases, occurring at J = J c ~ Tk, 
extends into the (B, J)-plane (Fig. 3(c)). The transition 
could thus be driven either by tuning J through J c at 
zero-field, or by tuning B through a critical B c at fixed 
J. With interlead charge transfer and L-R asymmetry 
both of which occur in our device, the transition line is of 
course broadened into a line of crossovers (Fig. 3(c)) with 
associated conductance peaks. For J J c in particular, 
where the system at zero-field is deep in the LS phase, 
B c ~ J ' I ' gun (and the Kondo screened phase for B > B c 
is naturally spin polarized, asymptotically approaching 
the free T_ state). This behavior is observed clearly in 
Fig. 3(d). In the center of the (1,1) regime, J ~ 0.34 meV 
(as estimated above) yields B c ~ J/gfiB = 2.9 T, in very 
good agreement with the measured conductance peak in 
Fig. 3(d) As detuning increases, J becomes larger as 
noted above, and B c is correspondingly seen to increase. 
The evolution of conductance as a function of magnetic 
field is in good agreement with the NRG calculations as 
shown in Fig. 3(b,d). The calculations also confirm (not 
shown) that the asymmetry in V sc \ observed (Fig. 3(a)) 
in the differential conductance with increasing detuning, 
results from the asymmetry Tl ^= Tr. For symmetric 
coupling strengths, or precisely at the center of the (1,1) 
charge regime where the system is electron-hole symmet- 
ric [e.g. the leftmost plot in Fig. 3(a)], no asymmetry is 
observed in the differential conductance. 

Finally, we show that by varying Vbar (see Fig. 1(a)) 
we can directly tune the coupling strengths between the 
quantum dots and their leads [32j. This is illustrated in 
Fig. 4(a) which shows that for decreasing T4 ar , the con- 
ductance at the center of the (1,1) charge regime strongly 
increases. At the same time, the double peak struc- 
ture, clearly observed for higher values of Vb ar , gradually 
merges into one broad peak. This behavior can be under- 
stood as an increase in the coupling strengths Fr,^, as 
illustrated by NRG calculations, see Fig. 4(b), and con- 
sistent with the increase in widths of the Coulomb block- 
ade peaks in the stability diagrams. On increasing r^ j?, 
and hence decreasing U/Tl^r, the device becomes less 
strongly correlated and charge fluctuations consequently 
more significant, thereby eroding and ultimately destroy- 
ing the double peak structure that is characteristic of the 
strongly correlated regime of the DQD. 

In conclusion, we have investigated a tunable two- 
impurity Kondo system in a strongly correlated carbon 
nanotube DQD, in which the full range of charge regimes 
is accessible, and which amounts to a clear experimen- 
tal realization of a two-impurity Anderson model. In 
the (1,1) regime, we have shown that the exchange and 
Kondo energy scales can be varied by tuning the dot en- 
ergy levels or by varying the dot-lead tunnel couplings, 
providing the possibility of observing the crossover be- 
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Figure 4: (color online) (a) Observed stability diagrams of 
the DQD for different dot-lead coupling strengths as varied 
by Vbar- (b) Left: differential conductance in the center of 
the (1,1) charge regime for different dot-lead couplings. The 
curves are offset in steps of 0.05 e 2 /h. Right: NRG calcula- 
tions, with the ratio Yr/Yl = 2 fixed and Yr as indicated 
(curves again offset in steps of 0.05 e 2 /h). 



tween local singlet and Kondo screened phases. While 
charge fluctuations in this device prevented observation 
of the crossover at zero field, we were able to observe it 
at finite field, indicated by enhanced zero-bias conduc- 
tance. The work is readily extended to carbon nano- 
tube DQDs coupled to superconducting [33| or ferro- 
magnetic [HI leads. This allows experimental access to 
rich phase behavior controlled by the interplay between 
Kondo, exchange, and superconducting correlations, and 
further highlights the potential of carbon nanotube quan- 
tum dots to investigate correlated electron physics. 
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The following appendices comprise the supplementary material to Ref. [l|. We provide the following: (i) details of 
device fabrication; (ii) additional experimental data on the dependence of the double quantum dot (DQD) conductance 
on the applied magnetic field, temperature and interdot and dot-lead tunnel coupling; and (iii) details of the model 
Hamiltonian used in the numerical renormalization group (NRG) calculations, and how the conductance is calculated 
from the NRG. 

DEVICE FABRICATION 

Carbon nanotubes were synthesized by chemical vapor deposition using a procedure similar to Ref. Q. To form 
the catalyst particles, ~ 0.5 /Ltg/cm 3 Fe(N03).9H20 (Sigma-Aldrich) was sonicated in isopropanol. The bare Si/Si02 
substrate (300 nm thermal oxide) was dipped into the solution immediately after sonication and dried by air blowing. 
It was then heated to 900 °C in a 1" tube furnace under a hydrogen gas flow of 400 seem, during which the iron nitrate 
was reduced to iron particles. Once at this temperature, 500 seem of methane (containing natural isotope ratios) was 
added to the flow for 8 minutes. After the methane flow was stopped, the substrate was cooled under the same 400 
seem hydrogen flow. Alignment marks were defined by electron-beam lithography, and scanning electron microscopy 
and atomic force microscopy were used to locate the nanotubes in relation to these. Further steps of electron-beam 
lithography and evaporation contact the nanotubes with 20 nm thick, 300 nm wide gold ohmic contacts, separated 
by 700 nm, which also form the outer barriers of the quantum dots. The final electron-beam lithography layer writes 
the 100 nm wide and 100 nm spaced gates, which are made from two 1.2 nm layers of aluminum, oxidized in air and 
capped with a further 20 nm of titanium and 6 nm of gold. The source lead is grounded by the virtual earth of a 
xlO 7 current to voltage preamplifier, and the drain has a small variable dc bias applied. The chip is mounted on the 
cold finger of a 60 mK dilution refrigerator and connected to the measurement apparatus with filtered low-frequency 
lines. 



MAGNETIC FIELD AND TEMPERATURE BEHAVIOR 

Figure SI (a) shows the evolution of the double quantum dot stability diagram with magnetic field in steps of 0.5 
T. We first focus on the widths of the hexagons (1,0) and (2,1). In the atomic limit (no coupling between the dots 
and their leads) the width of the hexagon increases proportional to gfisB for a magnetic field applied perpendicularly 
to the nanotube axis, where we assume 5 = 2 for nanotubes. In practice our device has a relatively strong coupling 
between the dot and leads, and a linear relation between the width of the hexagons and B is a good approximation 
only in the limit of large magnetic field. Bearing this in mind, we obtain the coupling strengths (or lever arms) 
between the gate electrodes and the quantum dots, and extract charging energies U ~ 2.5 meV for both quantum 
dots and an electrostatic coupling energy U' ~ 0.6 meV. 

The differential conductance as a function of magnetic field is shown in Figs. Sl(b) and (c) - which are an extension 
of Fig. 3 in the main text - for different positions in the stability diagram. Fig. Sl(b) shows the differential conductance 
as a function of magnetic field for seven positions in the (1,1) charge region as indicated by the circles in the B = 
diagram in panel (a). Fig. Sl(c) shows the differential conductance as a function of magnetic field in the center of 
the (2,1) charge region as indicated by the square in the B = diagram in panel (a). The maximum conductance in 
the center of the (1,1) charge region is observed at a magnetic field B ~ 2.8 T as seen both in Fig. Sl(a) and (b). 
The maximum conductance in the center of the (2,1) charge region is observed at B = T, as expected, with an 
approximately linear splitting of the Kondo resonance as the magnetic field is increased. 

The temperature dependences of the zero-bias conductance in the centers of the (1,1) and (2,1) charge regions 
are shown in Fig. S2. They can readily be explained with reference to Fig. 1(c) of the main text, noting that the 
dominant effect of increasing temperature is to broaden the low-energy conductance peaks. In the center of the (1,1) 
region (red points) there is a gradual increase of zero-bias conductance with temperature, corresponding to a filling-in 
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Figure S 1: (a) Double quantum dot stability diagrams measured at T ~ 60 mK. The magnetic field is increased from in 
steps of 0.5 T as indicated, (b) Differential conductance as a function of magnetic field at positions along the diagonal of the 
(1,1) charge region as indicated by the closed circles (•) in the B — diagram of panel (a), (c) Differential conductance as a 
function of magnetic field in the center of the (2,1) charge region as indicated by the closed square (■) in the B — diagram 
of panel (a). 
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Figure S 2: Measured temperature dependence of the conductance in the center of the (1,1) charge region (red points) and 
(2,1) charge region (blue points). 



of the conductance dip in the third panel of Fig. 1(c). In the (2,1) region (blue points) by contrast, the broadening of 
the single conductance peak (2nd panel of Fig. 1(c)) with temperature results in a gradual decrease in the zero-bias 
conductance. NRG results for similar parameters (not shown) display the same qualitative trends. 
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Figure S 3: (a) Stability diagrams of the double quantum dot for five different barrier gate voltages Vbar- (b) Differential 
conductance in the center of the (1,1) charge region (left), the (2,1) charge region (middle) and (2,0) charge region (right) for 
a series of barrier gate voltages as indicated. Curves are offset in steps of 0.05 e 2 /h in each panel, (c) Differential conductance 
as a function of magnetic field in the center of the (1,1) charge region for Vbar = —15 mV. 



VARYING DOT-LEAD AND DOT-DOT COUPLING STRENGTHS 



By varying the barrier gate voltage Vbar we were able to tune both the interdot coupling and the coupling between 
the dot and leads. Which of these two effects dominates varies between different electronic states, the microscopic 
details of which are not presently understood. Nevertheless, the effects are highly reproducible and can be used to 
investigate the competition between exchange and Kondo correlations on the quantum dots. 

Figure S3 shows the stability and differential conductance in the center of the (1,1) charge region as Vbar is varied 
- an extension of Fig. 4 in the main text. In this case, the main effect of varying Vbar is an increase in the coupling 
to the leads. This is deduced from the increase in width (FWHM) of the Coulomb blockade peaks in the stability 
diagram as well as the increased width of the Kondo resonance observed in the (2,1) charge region, see middle panel 
of Fig. S3(b). The double peaks in the center of the (1,1) charge region broaden with decreasing Vbar to form a single 
broad peak centered around V s d = 0, interpreted as a crossover to the weakly correlated regime where U/Tl,r ^> 1, 
see main text. 

A different dependence on Vbar is observed for the neighboring 'effective' (1,1) charge region, see Fig. S4. As 
expected, a double peak structure is observed in the differential conductance. However, in this case, the interdot 
coupling t is changed, resulting in a decrease of the exchange energy J with increasing Vbar and correspondingly 
smaller peak splitting, see Fig. S4(c). The change in t is also evident in the stability diagram, see Fig. S4(b), where 
it is observed as a decrease in the separation (and curvature) of the triple points. The smallest peak splitting in the 
differential conductance was observed for Vbar = 15 mV, see Fig. S4(c). No further decrease in t was observed by 
increasing Vbar beyond this value (not shown). 



MODEL HAMILTONIAN 



The canonical model used to describe tunnel-coupled DQDs is the two-impurity Anderson model [3j-|6|. We take 
this as the basic model for our NRG calculations, and supplement it with an interdot Coulomb repulsion/electrostatic 
coupling term (Htji), which inevitably exists, and which we find necessary to obtain good agreement between theory 
and the experimental stability diagram shown in Fig. 1(b). 

The Hamiltonian, H = + H-q + H\ ly b, consists of terms for the isolated leads and dots, and respectively, 
together with a hybridization term, -f/hyb, describing electron tunneling between dot-^ (y = L,R) and lead-i/. The 
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Figure S 4: (a) Stability diagram of the double quantum dot. Here the Si backgate voltage is set to Vt g = —1.8 V and Vbar — 
mV. The circle indicates the approximate center of an effective (1,1) charge region, different from that discussed in the main 
text, (b) Conductance around a triple point pair for the area indicated in panel (a) by the dashed lines for three different 
Vbar — 5, 10, and 15 mV. (c) Differential conductance in the center of the leftmost (1,1) charge region, i.e., position indicated 
by the circle in panel (a), for a series of Vbar- The increase of the peak splitting correlates with an increase of the interdot 
coupling t. 



leads, for which 

Hl = e k c lka c uk^ (!) 

consist of equivalent bands of non-interacting electrons with single-particle e nerg ies ek and spin a =t, I (where c^ kcr 
creates a cr-spin electron in state k of lead v) . As usual for metallic leads [7Hl0fl , we take the leads to be flat bands 
with half-bandwidth D. For the dots, 

Hv=^e„n vtr + tJ2(dl a d Ra + h.c.^ 

- * (2) 

+ U r M r M + U'n L n R , 

where n va = S va d ll<y is the number operator for cr-spin electrons on dot v, and h v = n u ^ + fi v ^. e„ is the single-particle 
level energy for dot-v, t the interdot tunnel coupling, U the on-dot charging energy /Coulomb repulsion, and U' the 
interdot Coulomb interaction. The hybridization term coupling the dots to the leads is given by 

ffhyb = E ^ + h - c ) . ( 3 ) 

with tunnel-coupling V v to lead-t/. The strength of hybridization to lead-^ is embodied in Y v = irV^p, where 
p = NJ {2D) is the (uniform) total density of states of the lead, and N the number of lead-orbitals (N — > oo 
in the continuum limit). The effect of a magnetic field is readily encompassed by adding Hb = —gpsBSz, with 
S z = i ^2 v (n^ — ^1/4.) the total z-component of spin on the dots, B the magnetic field strength, and g the electron 
(/-factor. We have taken g = 2 in comparing to experiment, but note that a slightly smaller value would lead to even 
better agreement between experiment and theory in Fig. 3 of the paper. 

In our calculations we use the parameters U, U', t and T u determined experimentally, and a temperature T ~ 100 
mK. Experimentally, the dot level energies are controlled by two gate voltages, V g L and V s r. One might naively 
anticipate e„ oc — V gv but, as is physically natural given the experimental set-up (Fig. 1(a)), there is some 'cross 
talk' between the side gate voltages, and to obtain good agreement between experiment and theory it is necessary to 
include some dependence of on V 9 r and of e_R on V 9 l- Specifically, we find 



e L oc -V gL - aV g R, e R oc -aV gL - V gR , 



(4) 
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with a = 0.4 giving the best agreement. 



NRG CALCULATION OF CONDUCTANCE 



We use the full density matrix (FDM) formulation (Til. [T3| of the NRG 043 (for a recent review see Ref . [Toh . keeping 
around 2000 states per iteration. To reduce discretization error, we average over different lead discretizations and 
calculate the dot Green functions via the self-energy method 13| . This enables us to work reliably with a relatively 
large discretization parameter of A = 9. 

To calculate the zero-bias (V sc i — 0) conductance, G c , of the model at finite temperature, we use an exact Kubo-type 
current-current correlation function approach, derived via linear-response theory [l4l [l5T | 



where (with ji, v E {L, R}) 



G c (V sd = 0,T) = — [a LL 



T^V^Vy lim 



Illl 



a RR — 



G M „(w) 



(5) 



(6) 



The correlation functions C m „(cj), which can be calculated directly via the FDM- NRG, are given by 



(lj) 



dwe~ iut (-i) 9(t) ( [j^t),M0)h) 



(7) 



where [a, b]+ denotes an anticommutator and 9(t) is the unit step function. The defined by 



3n 



(8) 



are current operators (modulo constants) , with / M oo- = 77^ X)k c Mkc the operator for the zero-orbital of the Wilson 
chain representation of lead fi — L, R 043- At T = 0, the zero-bias conductance eqsJSJIH] takes the simpler, physically 
intuitive form [l6| 



G c (V sd = 0,T = 0) 



-T l Tr\G lr {u = 0)| 5 



(9) 



FT i 

Here Glr{w) is the off-diagonal (or interdot) retarded Green function, Glr{u) > —i9(t){{d La .(t),d Ra ] + }, which 
may also be calculated via the FDM-NRG; and u> = denotes the zero-bias Fermi level. Results for the T = 
conductance obtained this way agree very well (as they should) with those calculated using eqs(S]|S]as T — > 0. 

At finite bias nothing exact can be said and, in practice, approximations must inevitably be made. Here we make the 
common approximation of bias-independent self-energies (which as judged by previous comparison to experiment 
[l9| appears quite successful), and approximate the finite- bias conductance by 



8e 2 r°° 
G c {V sd7 T) = —T L T R duj 

h J-00 



_ A %M_ {1 _ A) ?/-M 



duj 



du) 



\G LR (u},T)\ 



(10) 



where f v [u>) = 



_ \Ju-,i„)/k B T 



1] 1 is the Fermi function for lead v = L, R, with chemical potentials fi R = XeV s d and 



fil = —(1 — \)eV s d- The quantity A (g [0, 1]) is a measure of how the voltage bias is partitioned between the leads; 
A = 0.5 being perfectly symmetric, and A = or 1 meaning the bias is applied entirely on one lead. We take A = 0.8 
throughout, this value being found to give the best agreement with the experimental data (in particular for the slopes 
of the Coulomb blockade peaks, which are known to depend on A (20I 21 1). 
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